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6.3 



Abstract 

We consider a deformation of five-dimensional warped gravity with bulk 
and boundary mass terms to quadratic order in the action. We show that 
massless zero modes occur for special choices of the masses. The tensor zero 
mode is a smooth deformation of the Randall- Sundrum graviton wavefunction 
and can be localized anywhere in the bulk. There is also a vector zero mode 
with similar localization properties, which is decoupled from conserved sources 
at tree level. Interestingly, there are no scalar modes, and the model is ghost- 
free at the linearized level. When the tensor zero mode is localized near the IR 
brane, the dual interpretation is a composite graviton describing an emergent 
(induced) theory of gravity at the IR scale. In this case Newton's law of gravity 
changes to a new power law below the millimeter scale, with an exponent that 
can even be irrational. 
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1 Introduction 



It is a striking fact that four-dimensional (4D) gravity can be localized in AdS^ by 
tuning the bulk and brane cosmological constants pQ. It is even more remarkable 
that this five- dimensional (5D) model has a dual 4D interpretation via the AdS / CFT 
correspondence The gauge-gravity duality has made the warped gravity setup 

an attractive setting for studying aspects of strongly coupled gauge dynamics, from 
both the formal and phenomenological perspective. A feature particularly useful for 
low-energy model building is that the nonfactorizable geometry localizes not only 
gravity, but also fields of different spin, such as scalars and fermions |HJ E] • 

An important difference of scalar and fermion localization from gravity is that 
fermion and scalar zero modes can be localized anywhere in the 5D bulk. This 
unrestricted localization is achieved by introducing bulk and boundary masses, with 
the degree of localization directly depending on the bulk mass parameter |5] . This has 
been particularly fruitful for model building, where models of UV/IR brane-localized 
fields, corresponding to hybrid theories of elementary and composite particles in the 
4D dual, have opened new phenomenological possibilities for the supersymmetric 
standard model [TpJ [TT] and the Higgs sector [T3J EH- In fact, this unrestricted 
localization can also occur for gauge fields. While it is well known that gauge bosons 
are not localized in the bulk [TH ITS] , it is possible to tune bulk and boundary mass 
terms so that a U(l) zero mode can be localized anywhere [TB] . 

The graviton appears to be different, but it is natural to ask whether gravity can 
similarly be localized anywhere in the bulk. The AdS / CFT correspondence provides 
our primary motivation for studying the derealization of the graviton in the Randall- 
Sundrum (RS) scenario. If the graviton zero mode could be localized on the IR 
brane, this would suggest that in the dual theory the graviton is a composite CFT 
state whereby dynamical gravity only emerges in the infrared. The dual UV theory 
would then be a pure gauge theory with no propagating massless spin-2 particle. 
It has been suggested that in an "emerging" ("fat") gravity model, the graviton 
compositeness may alleviate the UV sensitivity of the cosmological constant [T7j . 
One may worry that derealization comes at the price of losing general covariance in 
the bulk, but as we will see, this does not preclude the existence of a massless tensor 
mode. Furthermore, the Weinberg- Witten theorem is avoided if gravity is induced at 
the quantum level [THj. Thus, there does not appear to be any obstacle to smoothly 
deforming the graviton wave function away from the UV brane, as can be done for 
other bulk fields, and one hopes that a similar mechanism can be implemented for 
gravity. 

In this paper, we study a modification of the linear Einstein equations and show 
that, indeed, the graviton can be localized on the IR brane. Much like for bulk 
scalar and fermion fields, our modification requires the introduction of bulk and 
boundary mass terms for the metric perturbations. However adding a graviton mass 
is technically more involved compared to the case of lower spin fields. For instance, 
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already in 4D a Fierz-Pauli ^Hj mass term for the graviton leads to the well known 
vDVZ discontinuity at the linearized level [20] (while other choices lead to ghosts). 
This is due to the fact that the longitudinal polarization of a massive graviton remains 
coupled to matter in the limit of zero graviton mass. In the warped 5D case, the 
extra polarizations of the graviton can be systematically studied by exploiting the 
symmetries of the 4D Minkowski slices. This leads to the definition of tensor, vector, 
and scalar fluctuations with respect to the 4D Poincare transformations. In this 
way, we show that by a special choice of the bulk and boundary mass terms there 
exist tensor and vector zero modes that can be localized anywhere in the bulk. A 
priori, one should also expect scalar modes but, remarkably, they are absent in our 
model. Thus, up to quadratic order in the action, we will see that our model is ghost- 
free. This is in contrast to many modifications of gravity where ghost perturbations 
typically arise in the scalar sector. In particular, a previous study of massive gravity 
in warped geometry, which only introduced boundary mass terms, concluded that 
ghosts were present in the spectrum [2*T] . 

Below, we present the structure and summarize the main results of this paper: 

• We begin, in Section 2, with a study of the bulk equations of motion for the 
graviton perturbation with a bulk mass term added. The solutions of the 
linearized equations of motion for the metric perturbations lead to massless 
and massive tensor, vector, and scalar modes. These are subject to boundary 
conditions on the two branes, which we discuss next. 

• Boundary mass terms are introduced and the junction conditions are derived, 
in Section 3, for all the metric perturbations. Some details are presented in 
the Appendices; in Appendix A we present the bulk gravity action to quadratic 
order. In particular, in Section 3 we also find that a massless tensor and vector 
mode is consistent with the boundary conditions and that the model is ghost- 
free. For any allowed value of the bulk mass parameter, there exist two possible 
choices of boundary terms admitting a zero mode (see Fig. 1). 1 The quadratic 
action of the graviton zero mode for the two branches is given in Appendix B. At 
the massive level, we also obtain the Kaluza-Klein spectrum. Most importantly, 
we note that the scalar spectrum vanishes identically. Details of the quadratic 
action for scalar perturbations and the relevant boundary conditions are given 
in Appendix C. 

1 Curiously, for zero bulk mass there are also two solutions admitting a zero mode. One of them 
is the original RS solution with zero boundary masses, while the other, requiring nonvanishing 
boundary mass terms, can be termed the specular RS solution: its zero mode wave function, ~ 
l/A 2 (z) and peaked on the IR brane, is the inverse of the original RS solution's wave function, 
A 2 (z); see Sections 3-5 and Fig.l [specular: (adj.) of, relating to, or having the qualities of a 
mirror] . 
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• We continue, in Section 4, with a discussion of the properties of the massless 
modes. We show that the tensor mode can be localized anywhere in the bulk 
and gives rise to a finite 4D Planck mass. This then motivates us to consider 
localizing the graviton on the IR brane and to interpret the force it mediates as 
an emergent infrared phenomenon since there is no propagating massless spin 
2 mode in the UV. In Sections 4.3, 4.4, and Appendix D, we study Newton's 
law at short distances and show that it exhibits a new power law dependence, 
that can even be irrational. The crossover to a V oc l/r^, (5 > 1 (and real) 
behavior is observable at distances shorter than the infrared scale — but much 
larger than the AdS curvature scale, unlike the RS model — which can be taken 
to be in the sub-millimeter range. Gravity at such distances is presently under 
experimental investigation [2*2]-|25|. 

• In Section 5, we discuss the 4D holographic interpretation of our gravity solu- 
tion. While it is necessarily speculative, even the remotest plausibility of a 4D 
dual picture inspires some confidence in the existence of a nonlinear extension 
of our linearized gravity analysis. 

— In the first of the two branches — the one smoothly connected to the RS 
model — the dual involves a strongly coupled 4D CFT coupled to a source 
graviton (Section 5.1). The coupling is relevant precisely when the zero 
mode is localized near the IR brane. As a result of this coupling, the energy 
momentum tensor of the CFT should somehow acquire a large anomalous 
dimension. The conformal symmetry is broken at the infrared scale where 
a large Newton constant is induced by the CFT. The massless zero mode is 
thus a mixture of the source graviton and a CFT composite. The dominant 
component of the zero mode graviton, in the case of relevant coupling, 
comes from the broken CFT. We show how the change in Newton's law 
can be explained in detail in the dual picture. 

— The holographic dual of the other branch of the gravity solution, described 
in Section 5.2, where the graviton is always localized near the IR brane, is 
even more mysterious. Holography suggests that it involves a CFT coupled 
to a massive source graviton (with mass of order the AdS curvature scale). 
Thus the observed massless graviton mode is essentially "emergent" at low 
energies due to the strong infrared dynamics of the broken CFT. 

• We conclude, in Section 6, with a discussion of the effects of nonlinear terms 
to our analysis whose details are beyond the scope of the present paper. 

Finally, let us comment on the relation of our model to the "fat" graviton idea of 
Sundrum [T7j, which was designed to turn off gravity at short distances (precisely 
at the scales close to the experimental limit [22]- [213) and alleviate the cosmologi- 
cal constant problem. In the picture presented here — which is the first quantitative 
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model of a "fat" graviton, where gravity in many cases is an "emergent" low-energy 
phenomenon — the gravitational interaction becomes stronger at short length scales. 
The stronger UV gravity may be a more general phenomenon, or be simply due to the 
particular setup: gravity in a slice of AdSs, being dual to a large- N CFT, has a tower 
of stable "mesons" (the graviton Kaluza-Klein modes) contributing to the gravita- 
tional potential at short length scales, and may just signify that the "fat" graviton, 
without the effects of the Kaulza-Klein modes, should be looked for elsewhere. 

2 Warped gravity with a bulk mass 

2.1 The Randall-Sundrum background solution 

In the absence of mass terms, warped gravity in five dimensions is governed by the 
following 5D action: 

S = J d 5 xy/^ (M 3 R -2A)-J2 J d * x v 73 ^ (M 3 [K] + A) . , (1) 

where M is the 5D Planck scale, A is a bulk cosmological constant, and the sum is 
over two boundary three-branes with brane tensions A,. The quantity [K] denotes 
the jump of the trace of the extrinsic curvature across the brane. 

The solution to Einstein's equations is a slice of AdSs where the fifth dimension is 
compactified on an orbifold S 1 /Z 2 of radius R with the Randall-Sundrum metric pQ: 

ds 2 = e-^ri^dxPdb? + dy 2 = A 2 {z){r] ilu dx lx dx v + dz 2 ) = gf B dx A dx B , (2) 

where < y < ttR is the "fundamental domain" (where the bulk integral (0) is 
computed), k is the AdS curvature scale, A = — 6k 2 M 3 , and the Minkowski metric 
7]^ has signature (— +++). The Latin indices (A, B, . . . ) label all the 5D coordinates, 
while Greek indices (fi, u, . . . ) are restricted to the 4D coordinates. We will work with 
conformal coordinates defined by z = (e ky — l)/k and A(z) = (1 + kz)~ l is the warp 
factor. At the orbifold fixed points z = and z\ = (e nhR — l)/k there are two 
three-branes, the UV and the IR brane, respectively. 

In the 5D bulk the action is invariant under general coordinate transformations. 
At linear order the infinitesimal coordinate transformations are: 

X M ^ X M + ^M {x) ; (3) 

where it is convenient to split the transformation parameters into £ M = (£ M + <9 M £, £ 5 ) 
with the condition = 0. In this way the coordinate transformations are separated 
into two scalar (£, £ 5 ) and one vector (£ M ) transformation parameters (with respect 
to the Poincare symmetry of the 4D Minkowski background). In a covariant theory 
these can be used to eliminate five degrees of freedom (d.o.f). 
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The metric perturbations, Hab around the background Randall-Sundrum metric 
9ab i correspond to fifteen degrees of freedom, and a useful way to parametrize them 
is: 

ds 2 = A 2 (z) [(1 + 2<p)dz 2 + 2 (B tli + B^dz dx^ 

+ ((1 + 2^)7]^ + 2E^ V + E M + V) dx^dx^ , (4) 

= (g { AB + h AB )dx A dx B , (5) 

where E^ v ) = d^E v + d u E fl . We see that the metric perturbations are divided into 
three sectors: scalar, vector, and tensor (with respect to the Poincare symmetry of 
the 4D Minkowski background). 

Specifically, the tensor mode h^ u , is taken to satisfy the transverse (d^h^ = 0) and 
traceless (h^ = 0) conditions. It is gauge invariant under the infinitesimal coordinate 
transformations 0, and being symmetric, it contains 10 — 5 = 5 d.o.f. 

The vector modes B^, and are both taken to be transverse (d^B^ = d^E^ = 0), 
and consequently contain 6 d.o.f. Under the coordinate transformation Q the two 
vector modes transform as: 

B,^B,-^, (6) 
En—tEp — ^, (7) 

where prime (') denotes differentiation with respect to z. Notice that there is a gauge 
invariant combination: 

B, ^B,-E'^. (8) 

In a covariant theory, the orthogonal combination can be eliminated with £ M using 
the coordinate transformation (|7jl. leaving only the three polarizations contained in 

K 

Lastly, the scalar modes ip, <p, B, and E, represent four real degrees of freedom. 
Under the remaining infinitesimal coordinate transformations they transform as: 

V>-^-^£ 5 , (9) 

0-0-£ 5 '-^£ 5 , (io) 

B^B-e-e, (ii) 

E^E-i. (12) 
There are now two gauge invariant combinations: 

^=^-^(B-E'), (13) 

^=^-^{B-E')-{B-E')' . (14) 
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These two modes ip, and represent two polarizations in the case of a covariant 
theory. 

For later convenience, we also define the gauge invariant brane positions. In 
general, a boundary brane can be at the perturbed position Zi + Q , and, as we shall 
see, the perturbation Q couples to the scalar perturbations of the metric. Under the 
infinitesimal change of coordinates (p5|). Q — » C«+£ 5 • Hence, we can form the invariant 
combination: 

C = C + B-E' . (15) 

In summary, the perturbations of the metric contain 5 + 6 + 4 = 15 degrees of 
freedom in the three sectors, corresponding to the degrees of freedom of a symmetric 
5x5 tensor. In the case of a covariant theory, such as the original Randall-Sundrum 
model, they are reduced to 5+3 + 2 = 10 polarizations. Not all of them are necessarily 
dynamical, since (as we will see) some polarizations vanish due to constraint equations 
that result from the equations of motion (equivalently, they can be seen as Lagrange 
multipliers in the original action). This parametrization of the perturbations will be 
also useful when we add mass terms to the action. 



2.2 Adding a bulk mass term 

Let us now consider adding a mass term for the perturbations Iiab = 9ab — S'ab > 
where g^ B is the background Randall-Sundrum metric. The 5D bulk action becomes: 



S = J d 5 x^ [M 3 R — 2A — M 3 k 2 g w MN g {0) AB {a h MA h NB + b h MN h AB )] , 

(16) 

where a and b are real parameters. The previous classification of the perturbations 
Uab is useful since modes belonging to different representations are not coupled to 
each other at the linearized level. This means that we can write the Einstein equations 
in the three sectors independently. For the covariant case, it is easiest to compute the 
equations in the form 5Gf} = 0, where 5G is the linear perturbation of the Einstein 
tensor. This equation holds as long as there is only a cosmological constant in the 
bulk (adding dynamical fields in the bulk gives rise to a nonvanishing STj^ 1 ) . 

If we now include the bulk mass term, then the action ()16|) leads to the bulk 
Einstein equations: 

5G B - + 2M 3 k 2 [ahi + bhS%] = , (17) 

where h B = g^ AC Ucb and h = h\. Let us now separately consider the nontrivial 
Einstein equations in each of the three sectors: 
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2.2.1 Tensor 

The equation of motion for the tensor modes h^ u is the transverse-traceless part of 
the \xv component of (|17p. and is given by: 



Dh 



h" 4- 3 — - W 



4ak 2 A 2 'h 



[IV 







where □ = —d 2 + d 2 . Notice that this equation does not depend on the b part 
of the mass term (fTBj) since the tensor mode is traceless. The solution of the 



equation of motion is obtained by a separation of variables h 



X,Z) 



where \3H llv {x) = rr^H^x), with m representing the mass of the four- dimensional 



Kaluza-Klein modes. The massless mode solution is: 



while the massive modes are: 



h$(x,z) = A- 2 (z) 



CiJ. 



2 Vl+a 



m r . 



kA(z) 



+ C 2 Yr 



2 I 2JT+E 



kA(z) 



(19) 



(20) 



where C\,C 2 are arbitrary constants. We will consider only values a > — 1 , which 
include the Randall-Sundrum case (a = 0), and, as we will see, provides a general 
and interesting phenomenology. In the limiting case (a = —1), the massless solutions 
are degenerate. Also note that in the limit a — ► these modes become: 



"fj.ii 



[C 1 + C 2 A(z)-*]Hj°)> 



X 



h^ n \x,z) 



A' 



C l J 2 



m r . 



kA(z) 



+ C 2 Y 2 



m r 



kA(z) 



(21) 
(22) 



which, together with appropriate boundary conditions (see below) smoothly repro- 
duce the Randall-Sundrum solution PQ. 

2.2.2 Vector 

The 5/i and \w components of ()17j) lead to the Einstein equations for the vector 
modes: 



UBn - 4 a k 2 A 2 



B[. + 3'-^B u + Aak 2 A 2 E u 



, 
. 



(23) 
(24) 



A' 

A^ 

These equations can be decoupled by eliminating £> M using (JBJ to obtain equations 
which depend only on B^ and E^. This gives rise to a second order equation solely 
in terms of B^. 



DB^ + B" - kAB ' - (4 a + 3) k 2 A 2 B : 



(25) 



where 124)) can be used to obtain from B^. As for the tensor mode we see that 
the vector modes do not depend on the b part of the mass term (jlfijl since the vector 
modes are transverse and do not contribute to the trace of the perturbation in the 
mass term. The equation (|25)) is again solved by separating the variables. When 
a ^ the massless mode solutions are: 



EW(x,z) 



1 



C 2 A{z 



,-l+2vl+o 



2ak L 



Cl (1 - v^T^) A(z)- 2 ( 1+VT ^) 



+c 2 (i + v / rT^jA(z)^ i -v^) 

while the massive mode solutions are given by: 



A' 



2 y/l+a 



kA{z) 



+ C 2 Y, 



2 Vl+a 



kA(z) 



b ( ; ] (x) 



(26) 



(27) 



(28) 



where C\,Ci are arbitrary constants and E^ are obtained from (J24*|) . 

We can compare this with the massless case. When a = the equations of motion 



and 



only depend on From 



we see that the vector mode is always 



massless and since 
given by: 



is a first order differential equation there is only one solution 



Bf(x,z)=C 1 A-\z) bJx) 



(29) 



where C\ is an arbitrary constant. We will see later that the boundary conditions 
will eliminate this mode, as is well known for the RS case. 



2.2.3 Scalar 

The scalar equations are obtained from the 55 , 5/j, and \iv components of (|17j). This 
leads to the following equations involving the scalar modes: 

4 \ V ~ I^j + I ^ \ a< P + b{Ai;+ <j) + DE)] 

v A 3 

+ 2^ - 4a k 2 A 2 E 

+ ^ k 2 A 2 [a^ + &(4^ + + UE)\ 

Notice that unlike the vector and tensor modes the coefficient b now appears in the 
scalar equations. This is because the scalar modes do contribute to the trace of the 
metric perturbations. 



= o, 


(30) 


= o, 


(31) 


= o, 


(32) 


= . 


(33) 
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When a = b = 0, only the gauge invariant combinations of the metric perturba- 
tions appear in the Einstein equations. The system of equations is straightforward 
to solve and from ()31|) and (|32|) we obtain: 



= -2 if> 



i>{x,z) = C x A- 2 (z)S{x) 



(34) 



where C\ is an arbitrary constant. Hence, as in the vector sector, only one mode 
is present in the RS case. Using equation (|3T?|) this mode (the radion) is massless, 
\3S(x) = 0. The remaining equation ()33|) is degenerate and consistent with the 
solution flHU). 

The solution of the system of equations (|3lH) - (j3^j) for the scalar perturbations 
ip,4>,B and E when a / and b ^ is more involved. As first step we will eliminate 
ip and in terms of ij) and 0. Thus, using flTB^l . (fPfj) . (pITj) . and (J3~2|) we obtain: 



E 
B 



Aak 2 A 2 
3 

2ak 2 A 2 

9- 1 



AakA 
1 



Aak 2 A 2 



+ 2^ 
i>' + kA <p 

-ft + 4kA [4>- r 
Ai)" + 3kAft 



(35) 
(36) 
(37) 
(38) 



Using these equations the remaining two scalar equations, (J30|) and (j33|l can then 
be written in terms ip and 0. This leads to two coupled second order differential 
equations which can be solved for tp and 0. These equations are equivalent to a 
fourth order differential equation, which is difficult to solve in general. However, the 
coupled differential equations magically simplify for the "Fierz-Pauli" choice: 2 



b = —a 



(39) 



of the bulk mass parameters in (JTHJ), since the sum of the two coupled equations only 
involves the combination ip — (p. The orthogonal combination ip + is then most 
easily obtained from ()30|) . Therefore, defining the linear combinations: 



the remaining two scalar equations are: 



Y = ?/> + 



□ X + X" + kA X' - 4 (1 + a) k 2 A 2 X 
1 

Y 



(3 + 4 a) k 2 A 2 



2kAX' - (5 + 4 a) k 2 A 2 X + -DX 



(40) 

(41) 
(42) 



2 In fact the major motivation for this choice is provided by the presence of (DE) 2 in the quadratic 
action for the scalar perturbations. As in the 4D case these terms cancel for the Fierz-Pauli 
choice; see Appendix C. 
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Hence, we have obtained a second order equation in terms of X only. The remaining 
equation is an algebraic equation for Y, which is trivially solved in terms of X. It is 
straightforward to obtain the general solution of the remaining differential equation 
(JJTJ). Again separating the variables and writing X(x,z) = f(z)S(x) we find that 
the massless modes solving these equations are: 



X(°\x,z) 
Y<P\x,z) 



C x A(z) 



-2Vl+a 



C 2 A(z 



-1 



4a + 3 



d 2V1 + a- 1 A(z 



x 



-2-JT+a 



+ C 2 foVT+a + lY 'A(z, 



2Vl+a 



(43) 



S (0 \x) , (44) 



where Ci, C% are arbitrary constants and the four dimensional mode obeys □S ,< -°- ) (x) 
0. If we now substitute these general solutions back into (f3l) |) -()38 |) we obtain: 



E(°\x,z) 



B^(x,z) 
^\x,z) 
<P^\x,z) 



2a (3 + 4a) k 2 



Ci(2 a + 3(1 - VTTaj) A(z)-^ 1 -" /TTE) 



+ C 2 (2 a + 3 (1 + vTTa)) 4(z) 



-2(l+v / T+^) 









(45) 

(46) 
(47) 
(48) 



Remarkably, the bulk equations of motion have forced all the massless scalar per- 
turbations to become zero, except for the mode! Furthermore, as shown in 
Appendix O this mode gives a vanishing contribution to the action for the perturba- 
tions, and therefore is not physical (at least, to quadratic order in the perturbations). 
Thus, the addition of the bulk mass term with the Fierz-Pauli choice ()39)1 has com- 
pletely eliminated the massless (radion) mode. 
The massive mode solutions are: 



X {n \x,z) 



2 Vl+a 



kA(z) 



&2 ^2vTP^ 



m r . 



kA(z) 



S (n) { 



X) 



(49) 



where US {n \ 



x 



m 2 S^(x) 



and (x,z) is obtained from (J42|) . As we shall see, 



also these modes disappear, once the boundary conditions are taken into account. 



3 Brane localized mass terms 
3.1 Boundary conditions 

In order to determine the mass spectrum from the general solutions we need to 
derive the boundary conditions satisfied by the bulk modes on the branes. We will 
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also add boundary mass terms on the branes since this will be crucial for obtaining 
a deformation of the RS solution. Hence, consider the following brane action at the 
location zc 



ASi = -kM 3 J d 4 x V-70 h a(3 (a { 7^ 7^ + $ 7^ ^ 



(50) 



where 70,^ = A 2 rj^ is the background induced metric on the boundary, with A 
evaluated at the (unperturbed) location of the brane, and are the perturbations 
of the induced metric, 

liiu = lQ,nv + . (51) 

The total induced metric 7 is evaluated at the perturbed brane position Zi + Q (x^) . 
The brane displacements Q constitute two additional scalar modes of the system 
which have support only on the two boundaries. It is worth noting that the action (|5U|) 
is not the unique boundary action which one could choose. This term itself is not 
covariant, and there are several inequivalent possibilities (as we discuss below in 
Section EPj) . Our choice was motivated by the fact that (|5U|) is rather simple and 
natural. 

The boundary mass terms give a contribution to the energy-momentum tensor on 
the boundary: 

5S liP = -^= 6 -^ = -AkM\ah^ + (3h l0llv ) , (52) 



where to linear order in the perturbations: 



hni, — lA 2 {zi 



( A'(z-) \ 1 / ~ 



(53) 



and h = M = ^ h vlx , with 7^ the inverse background induced metric. This 
contribution must be added to the standard RS piece arising from the brane tensions, 
namely Sj$ = — Xi^^u, where Aj = ±6kM 3 and ± refers to the UV/IR brane, 
respectively. Hence the total energy-momentum tensor is: 

S,u = S$ + 6S, U . (54) 

This expression is related to the jump of the extrinsic curvature across each 
brane. Note that the extrinsic curvature is computed from bulk quantities, and 
formally it is not affected by the bulk mass term in (|16p. The extrinsic curvature, up 
to first order in the perturbations, is given by: 

= V^n u = d^n u - T%, n 5 = d^n u - 1%,A (l + + ^(^j . (55) 
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Evaluating it from the bulk geometry gives: 



K 



5K, 



A' r]^ + 5K^ , 
2 A' 




, A' ± A" , 

ib 6 H C 

A A A 



E' + —E-B 
A 




The junction conditions are then: 

M 3 



-5. 



\w 1 



(56) 



(57) 



where K, 



K, 



. } means the jump across the brane (with Z 2 
symmetry imposed). Evaluating the junction conditions explicitly in the three sectors 



Kh^y and 



gives: 




tensor 


K- 


vector 




scalar 


: il/ - 


scalar 


E' — 



= ±4aikA , 

T^aikAE^ , 
A' 
~A 
B- 



= ±AkA 
d = ±4(XikAE 



A' \ ai + Pi 



+ UE 



(58) 
(59) 

(60) 

(61) 



These expressions are to be evaluated at the brane locations z (upper sign) and z\ 
(lower sign). They are valid in general, both for the massless and massive modes, 
and we will use them to determine the mass spectrum. 

In the scalar sector, there is an additional boundary condition, which is obtained 
directly from varying the quadratic action for the scalar perturbations with respect 
to the brane displacement Q (see Appendix O for details). When the mass terms are 
absent, this equation is actually redundant with respect to the boundary conditions 
given above. This is not surprising since it is a consequence of general covariance 
of the model. In the present case, the inclusion of the mass terms removes this 
degeneracy, leading to the additional boundary condition. For b = —a ^ , and 
fa = —OLi 7^ , the additional boundary condition can be rewritten as: 

scalar : 4 ip + UE = , (62) 

which, like the previous boundary conditions, only holds at the two boundaries. 

3.2 Tensor modes 

Applying the boundary condition (|58J) to the tensor mode general solution (|T9|) gives: 



(1 - v / TT^ =F 2on)Gx A( Zi ) 2VTT " + (1 + VTT^ =F 2a,i)C 2 A( Zi y 2VT ^ = 



(63) 
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Figure 1: The range of bulk (a) and brane (a) mass parameters that lead to a 4D 
massless graviton. There are two branches, oe±, joined together at the limiting value 
a = — 1. The RS model (a = a = 0) is a special case on the a_ branch. Its mirror 
image on the a + branch is the "specular" RS solution (a = 0, a = 1). 

where this equation is evaluated at z (upper sign) or z 1 (lower sign). For generic 
mass parameters a, on the only solution is C\ — Ci — 0, and there is no massless 
graviton. However, when: 3 

ao = — cxi = a , a = -(1 =F VT+~a) = a T , (64) 

the coefficient C\ [C2) drops from the boundary condition, and (jfiHj) simply gives 
C2 = (Ci = 0). Hence, the massless tensor mode (fTHj) becomes: 

h fll/ (x,z) = N T A(z)- ia H^(x) , (65) 

where Nt is the overall normalization constant which is determined from the quadratic 
action in the perturbations and cannot be determined by the boundary conditions. 

This form of the solution is only meaningful when the bulk mass parameter a > 
— 1. This corresponds to a = a + (aJ) for a > 1/2 (a < 1/2), so that the full range of 
the boundary mass parameter a is covered by the two branches a±. This behavior 
is plotted in Figure [TJ When a, a_ — > 0, the bulk and boundary mass terms become 
zero, and the a_ mode reduces to the usual RS tensor mode that is constant in the 
z coordinate. Hence the a_ mode is a smooth deformation of the RS tensor mode 
from a = to —00 < a < 1/2. On the other hand the a + mode can only exist when 
the boundary mass is nonzero and corresponds to the deformation of the RS tensor 
mode to values 1/2 < a < 00. Later, we will see that this corresponds to localizing 
modes continuously from the UV brane to the IR brane. 

3 Since we restrict ourselves to the Fierz-Pauli choice $ = — a, , we also set (3q — ~Pi = P f° r 
the scalar modes. 
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It is worth highlighting the presence of the massless zero mode for a = , a + = 1 , 
arising for a special value of the brane mass parameter, and without a mass term 
in the bulk. In this case, general covariance is broken only at the location of the 
two branes, while the 5D bulk is general coordinate invariant. As is clear in Fig. Q 
this model is the "specular" version of the RS mode. Indeed, as we have seen, the 
bulk equation is a second order differential equation, and in particular the value 
a = leads to two linearly independent bulk solutions. The choice of a + = 1 leads 
to the survival of the massless mode that is normally removed by the RS boundary 
conditions (and vice versa). 

For the massive solutions (J2*U|) the boundary condition (J55J) for a = a± can be 
written as 

J , ( _jnn_\ J 



2Vl+^±l \kA(z Q ) J °2VT+E±1 I kA(z 1 ) 



J- 9. /TXTT-i-1 i. a i - \ I J- 1 



(66) 

2 v / l+a±l ykA(zo) J 1 2y/l+a±l ykA(zi) ) 

Since A(zq) = 1 and A(zx) = e~ nkR we obtain in the limit ke~ nkR m n k the 
Kaluza-Klein mass spectrum 

m n ~ (^n + y/TT^ - |J vr ke~* kR , n=l,2,3,.... (67) 

This approximation for the mass spectrum becomes increasingly better as n grows. 
When a = we recover the RS Kaluza-Klein mass spectrum pQ. 



3.3 Vector modes 

Since the massless tensor mode only exists for a = a± we will also assume this value 
of a for the vector mode. The boundary condition (jHnj) for the massless vector mode 
(f2T)j) becomes: 

d(l + -a a.) A( Zi )-( 1+2 ^ + C 2 (l + -a a+) A^ 1 - 2 ^^ = , (68) 
a a 

where this condition is imposed at Zq and z\. When a = a_, the term proportional 
to C*2 in ()68|) vanishes (since = —a/4), which then simply enforces C\ = . 

On the contrary, for a = a + the boundary condition (j^Hj) leads to C2 = 0. Thus the 
massless vector mode solution becomes: 

B®(x,z) = N v A{z) l -' a bf{x) , (69) 

Ejp{z,z) = ^(l-a)A(z)-^b^(x), (70) 

where Ny is a normalization constant and a = a± for a > 1/2 (a < 1/2). Notice 
that when a, a_ — > the mode which gets killed, C\ = , is precisely the solution 
which would have smoothly gone to the RS vector mode (|29p. Hence, for a < 1/2 
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(including, in particular, the RS limit a = 0) this mode is always killed by the Z 2 
symmetry. However, we see that for nonzero boundary masses a^O there is a new 
massless mode present, given by (JHHJ), that is instead absent in the RS case. Just like 
the tensor mode, this mode can be localized anywhere in the bulk. 

The boundary condition (|59p for the massive solutions ()28|) leads to the same 
equation as for the graviton tensor modes, and hence to an identical Kaluza- 
Klein mass spectrum (J67|) . This is not too surprising since both modes originate from 
the five-dimensional graviton. 

3.4 Scalar modes 

Let us discuss the zero and massive scalar modes separately. For the zero modes, 
we have seen that the bulk equations enforce i/j = = B = . The boundary 
condition (|6L)jl then leads to Q = for both branes. This leaves only the scalar mode 
E which can in principle be nonvanishing. However, an explicit calculation shows 
that if E is the only mode present, it does not contribute to the quadratic action for 
the perturbations. Hence we conclude that, at least at the linearized level, there are 
no massless scalar modes in the theory. 4 

For the massive modes, there are three nontrivial boundary conditions (J6Uj) . (J61|) . 
and ()62|) . Since each equation holds on both branes, we have a system of six equations 
in four variables for each massive mode. The variables are the displacements of the 
two branes Co,i and from the mode the ratio C\jCi and the mass m n . Since 
this is an overdetermined system of equations we can only hope to have nontrivial 
massive scalar modes if some of these equations are degenerate. 

To check if there is a degeneracy first note that, after some algebra, Eq. (|62|) forces 
the mode in to vanish at both branes. The remaining equations can then 
be combined to eliminate Q, and together with (p)2"j) . they also force X' to vanish at 
the two branes. Thus there is no degeneracy and the requirements of: 

X (n) = X (n)' _ q ^ at both boundaries ; ( 71 ) 

can only be satisfied for C\ — Ci — . This immediately implies that the massive 
scalar modes are also absent. 

The absence of massive scalar modes is due to a "mismatch" between the bulk 
and boundary mass terms that we have introduced. Both terms can be considered 
as massive deformations of the original Randall-Sundrum proposal. We have previ- 
ously shown that the deformation (|52j) of the boundary action "matches" with the 
deformation (J 16)) of the bulk action to produce nontrivial tensor and vector modes. 
However, we now see that this is not the case for the scalar sector. It is possible that 
some inequivalent choice for the boundary actions can also accommodate nonvanish- 
ing scalar modes. As we already mentioned, the action (p)2"|) is not covariant, so other 

4 We have verified that for the "specular" solution, a = 0, a = 1, where the calculation is slightly 
different, there is also no scalar mode in the spectrum. 
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inequivalent actions can be considered. This is particularly relevant for the scalar 
modes, which have the additional ambiguity (with respect to the other sectors) of 
the brane positions. For instance, the perturbation of the induced metric in (|52|) is 
defined as: 



where 7^ is the (total) induced metric, with scalar perturbations included. In this 
way, the background induced metric is evaluated at the unperturbed brane positions 
£0,1 ■ At linear order, this leads to the expression (|53)1. Alternatively, one can choose 
to evaluate also the second term in (fT2~j) at the perturbed brane positions Z{ + Q. 
This results in omitting the terms proportional to Q in the expression ([53)1 . As a 
consequence, eq. (160)1 also appears without the two terms proportional to Q, while 
clearly the vector and tensor sectors are instead unaffected by this choice. We have 
verified that this choice also leads to vanishing scalar modes. Another possibility 
would be not to perturb the brane positions at all. This then leaves two fewer 
degrees of freedom, together with two fewer equations (i.e. Eq. (|B^)l evaluated at the 
two branes). Again we have verified that in this case the remaining equations are not 
degenerate, so that no scalar modes are present. 

In summary, all the inequivalent choices for the brane mass terms that we have 
considered lead to vanishing scalar modes. Although these choices appear as the most 
natural possibilities, clearly we have not exhausted all the possible choices. Hence, we 
cannot rule out the possibility that some other brane term can lead to nonvanishing 
modes also in the scalar sector. 

4 Mode properties 

4.1 Localization of the zero modes 

The existence of zero modes follows from the relation (|64|1 between the bulk and 
boundary mass parameters. However, this relation only fixes one of the parameters, 
say a, so that there is still freedom to choose a. Since the wavefunction depends on 
a we can arbitrarily localize the zero modes anywhere in the bulk. This is similar 
to previous studies [§] involving fields with spin < 2. Consider the effective four- 
dimensional action for the tensor modes. In terms of the RS variables y we obtain: 



where we have used the solution (|63jl in the effective bulk action at quadratic order 
(see Appendix |B] note that indices in (|75jl and ((75)1 are raised with 77^). From ((75)1 
we see that the tensor zero mode wavefunction has a y dependence: 




(72) 



S eff = ——MNr 




Jo 



fr iv) oc e 



(1— Aa)ky 



(74) 
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When a = we obtain the RS tensor mode localized on the UV brane. However we 
now see that by varying a we can smoothly deform the tensor mode to be localized 
anywhere. For a < the mode becomes even more localized on the UV brane, 
compared to the original RS scenario. On the other hand, for a > the tensor mode 
is delocalized away from the UV brane towards the IR brane. The transition occurs 
for a = 1/4 where the tensor mode is completely flat. As a becomes larger than 
1/4 the tensor mode becomes more and more localized on the IR brane. Hence, we 
have a continuous deformation of the original RS tensor mode from being completely 
localized on the UV brane to being completely localized on the IR brane! 
Similarly for the vector mode we obtain the effective action: 

S eff = ~M 3 N^ J d 4 x £ R dye~ 4 ^~ 2a)ky F$(x)F® fU '(x) + . . . , (75) 

where we have used in the quadratic bulk action and F$/(x) = d,j0\x) — 
d v bp{x). Therefore, the vector zero mode has the wavefunction dependence: 

/f^oce- 2 ^ 2 ^. (76) 

Again we see that the vector mode can be localized anywhere in the bulk. In partic- 
ular when a = 1/2 the vector mode is flat. Note that there is a difference between 
the tensor and vector mode wavefunction dependence so that these modes are never 
localized at the same place. 



4.2 4D Planck mass 

The 4D (reduced) Planck mass (Mp) is defined by assuming that there is matter on 
the branes described by an energy- momentum tensor T$\ 

Smatter = J d*X N T T® (x) , 

= Jd 4 x^-T^^H^(x), (77) 

where H$)(x) is the canonically normalized tensor zero mode. Thus, the four dimen- 
sional Planck mass is: 

= M^ A^( Zi ) r 2{Aa _ l)nkR _ ^ 
k (4a - 1) L J 

When a = this expression reduces to the RS result [T]: 

Ml ^ \ , (79) 
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assuming irkR 3> 1. In particular notice that the RS result does not depend on where 
the matter is localized. 

However, when a^O the Planck mass depends on whether matter is localized on 
the UV brane or IR brane. For matter on the UV brane, A(zq) = 1, and we obtain: 

J M 3 2(4a-l)wkR n v. 1 

fe(4a-l) e " ^ 4 ' 

M 3 2ttR a = \, (80) 



k{Aa - 1) 



M3 a < 1 



fc(l-4a) 4 



where we have assumed irkR 3> 1 in the approximate expressions. These expressions 
are consistent with the behavior of the graviton wavefunction. In particular when 
a < 1/4 the Planck mass is consistent with the fact that the graviton wavefunction 
is localized on the UV brane and there is no exponential expression. The a = 1/4 
result follows from the fact that in this limit the graviton wavefunction is flat and we 
recover the flat space result |2Zj. For a > 1/4 the graviton is localized towards the 
IR brane and the coupling to matter in the UV brane is exponentially suppressed. 
Instead when matter is localized on the IR brane, A(z\) = e~' KkR ) we obtain: 

e -2nkR a > 1 



M 3 fc(4a-l) " " " 4 ' 

2 _ 1V± \ a -2irkR -8airkR-\ ) ]\/f3 o*n- T? ^ — 2nkR „, _ 1 



Ml = [e- 2 * kR - e - 8a * kR ] ~ <^ M 3 2irRe-^ kR 



4 ' 



M 3 e - 8a * k R a K 1 



fe(l-4o) ° " ^ 4 ' 



again assuming nkR 3> 1 in the approximate expressions. These expressions are 
consistent with the localization properties of the graviton, except that now there is 
a nontrivial warp factor for the matter on the IR brane. 

An interesting phenomenological scenario occurs when the graviton is localized 
on the IR brane and all the standard model matter is located on the UV brane. If 
we associate the IR brane with the millimeter scale 10~ 3 eV, then for M ~ k ~ TeV, 
we obtain the usual Planck mass Mp for a = 1/2. The scale of the UV brane is the 
TeV scale. In the bulk the tensor mode has a profile e hy , and it is therefore localized 
away from the UV brane. This explains the weakness of gravity with respect to the 
gauge interactions. One can also verify that the same exponential suppression also 
holds for the coupling of the massive tensor modes. The vector zero mode is instead 
flat in the bulk. However, the vector mode does not couple to brane sources with a 
conserved energy-momentum tensor. Hence, conserved matter on the brane interacts 
only with the graviton tensor mode. 

We will see that the dual picture is particularly interesting since it indicates that 
the graviton is composite at the millimeter scale, and emerges from the dual CFT as 
a massless bound state. This emergent gravity picture is not in contradiction with 
the Weinberg- Witten theorem ^H] because as in models of induced gravity, there is 
a source gravitational background which invalidates the assumptions of the theorem. 
The corresponding energy-momentum tensor is no longer conserved and can obtain 
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a (large) anomalous dimension. Since the UV scale is associated with the TeV scale 
the gauge hierarchy problem is trivially solved (as usual, additional dynamics may 
be necessary to generate and stabilize the IR scale). 



4.3 Short range modifications of gravity 



We are interested in the gravitational interaction between nonrelativistic matter 
sources on the UV brane. This interaction is due to the tensor perturbation of the 
geometry. It includes the contributions from both the zero mode (already studied in 
Section 4.2 above) and the massive Kaluza-Klein modes. As we have seen in ([67)1 . 
the lowest Kaluza-Klein mass is at the IR scale mi ~ k A(z{) . This scale sets the 
distance above which gravity is standard. We are interested in the gravitational in- 
teractions at distances r smaller than this, but still much greater than 1/k (so that 
they can be phenomenologically relevant): 



1 1 

7- < r < — - 
k kAi 



(82) 



In 



where we have denoted A\ = A(z\) . The computation is given in Appendix D 
the following we present the final result and discuss it. 

The contribution from the Kaluza-Klein modes presents two distinct behaviors 
depending on whether a is smaller or greater than 1/4 . For either region, we define a 
positive parameter £ = |4a — 1| (therefore, the two regions are joined at £ = ; notice 
however that we have chosen £ to be positive in both of them). The gravitational 
potential is found to be: 



V(r) 



V{r) 



fi 



Mpr 



1 + 



2 T (2f) 1 
WW) ~{2kT^ 



£ = 1 - 4a , a < 1/4 , 



fi 



1 + 



2T(2g) 1 

(2fcA 1 r) 2 « 



M 2 p r 



x 



2T(2g) 1 

T^m (2kA ir ) 2 t 



r £ 



r ^ 



;f = 4a -l,a> 1/4.(84) 



l 



The 1/r term in these two expressions is the contribution from the zero mode, which 
reproduces the standard Newtonian gravity at large distances. Instead, the second 
term represents the interaction mediated by the Kaluza-Klein massive modes in each 
case. 

Note that there is a striking phenomenological difference between the two cases, 
shown in (JHHj) and (|84p. respectively. Indeed, in the first case a < 1/4, the corrections 
to the Newtonian potential are relevant only at distances r ^ k~ l , that is only at 
the UV scale. 5 However, in the second case of a > 1/4, the gravitational potential 

5 The RS model, characterized by a — belongs to this interval; a — corresponds to £ = 1 , so 
that we recover the 1 /r 3 correction computed in [2Hj • 
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is strongly modified already at the much larger (and, possibly, phenomenologically 
relevant) IR scale, r ~ (A\k) 1 . Moreover, below this distance the 1/r term dis- 
appears and the potential is given solely by a power law l/r 2 ^ +1 . The cancellation 
of the 1/r term will be explained in the next section, but it is already clear that 
the corresponding force is stronger than the usual gravity. So, for a > 1/4 standard 
gravity only emerges at infrared distance scales r ^ (Ax k)~ . 

Actually, this behavior is due to the different localization of the zero mode in 
the two regimes. 6 For a < 1/4, the zero mode is localized towards the UV brane, 
and the relative contribution from the Kaluza-Klein modes can be neglected. On 
the contrary, for a > 1/4 the zero mode is localized towards the IR brane and 
the relative contribution of the massive modes significantly increases becoming the 
dominant contribution for r ^ (Ax k)^ 1 — the regime indicated in (|82j) . We will see 
that for a > 1/4 this behavior is consistent with the graviton being a composite at 
the IR scale. 

Let us briefly examine the experimental consequences of our composite graviton 
model in the intermediate regime ()82|) . Assuming that the IR (or compositeness) 
scale is related to the cosmological constant then the IR scale is ~ 1CT 3 eV. As we 
have seen above at energies smaller than this scale we have the usual Newton law. 
The striking experimental signal of our model would be that Newton's law of gravity 
changes to a new power law r _1 — > r~ 2 ^~ l below ~ 0.1 mm, which could be fractional 
or even irrational! 



4.4 The composite graviton and Green's function analysis 



The modifications of Newton's law and the compositeness of the graviton can also be 
studied by computing the 5D bulk Green's function associated with the tensor mode 
h^. The formalism for computing 5D propagators in a slice of AdS can be found 
in Appendix A of Ref. [TU]. The Lorentz index structure of h^ v can be neglected, 
and then the Green's function of the tensor mode is identical to that of a bulk scalar 
mode with mass M| = 4ak 2 and boundary condition r = 4a± (using the notation 
of Ref. JUJ). A straightforward substitution into the general expressions in Ref. [TU] 
then gives rise to the Planck brane-Planck brane Green's function: 



G(x, z ; x', z Q ) 



d 4 p 



where: 



Gp( z o, Zq) 



Ivk 



1 

P 



I v (pL)K v±1 (p/k) + K u (pL)I u±1 (p/k) 



55) 



6) 



h(pL)K u (p/k) - K v (pL)I v (p/k) 

6 We do not expect the potential to be discontinuous at a — 1/4. However, the two expres- 
sions H8H(I and 1)84(1 become invalid for £ — > 0, since the approximations adopted to derive them 
break down in this limit. We expect that in this region the exact potential quickly (but smoothly) 
interpolates between the values given on the two sides. 
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Here p = y/p*, v = ±(4a± — 1) = v±, L^ 1 = kA(zi) = ke~ nkR is the IR scale, and 
I„(z),K v (z) are the modified Bessel functions. The Green's function f!86j) has been 
written as a sum of two terms [2H] • The first term represents a part of the zero mode 
contribution, while the more complicated second term receives contributions from 
both the zero mode and the Kaluza-Klein tower of massive modes: the zeros of the 
denominator of the second term at p — im n precisely reproduce the Kaluza-Klein 
mass spectrum (|fi7j) . Moreover, the contribution from the zero mode is clear from 
the fact that also this term scales as 1/p 2 for p — > 0. 

The contribution of the massive states can be simplified by expanding the Green's 
function in the 4D momentum p. Since k is the AdS curvature scale we will assume 
p <C k. Consider first the case of a < 1/4. This corresponds to the localization of 
the graviton on the UV brane. In the limit of either pL ^> 1 or pL <C 1 the dominant 
contribution to the Green's function is always: 

k 2M 3 1 

G p (z , z ) ~ 2(4« _!)_ + ...- -__ + . . . , (87) 

where in the second equality we have substituted for the Planck mass from (|8U|). 
Thus, we see that at low energies pL <C 1 only the zero mode contributes to the 
Green's function and the massive modes are decoupled. Similarly, for pL ^> 1 the 
zero mode is again dominant. This behavior is consistent with the gravitational 
potential (|83|). and the fact that the tensor zero mode is localized on the UV brane 
and remains pointlike for all energies p <^ik. 

Let us now consider the case where the graviton is localized on the IR brane. 
Assuming a > 1/4, we find that in the limit pL <C 1: 

G p (z , z ) ~ -2ie~^ kR - 2 + ■ ■ ■ ~ 1 + . . . , (88) 

p 2 Mp p 2 

where £ = 4a — 1 and the Planck mass has been substituted from (|SUJ). Again we see 
that the zero mode gives the dominant contribution at energies below the IR scale. 
There is now also an extra volume suppression that is absorbed into the gravitational 
coupling. On the other hand at energies above the IR scale, when pL ^> 1, we 
obtain: 7 

G r Uz 0) K (89) 
p K^p/k) 



1 



2(1-0* 



p 2 _ (p_\^-i r(2-p 

4A*(£-l)(£-2) \2k) r(0 



Remarkably, there is now no longer a dominant zero mode contribution since the 
massive Kaluza-Klein modes cancel the leading 1/p 2 contribution in ([86))! We will 

7 This expansion assumes £ ^ 1/2 . For £ = 1/2 , the analytic terms are not present in the series 
expansion of 
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see that this behavior is consistent with the fact that the graviton is composite at 
the IR scale L -1 . When pL <C 1 the graviton is essentially pointlike and contributes 
in the usual way, but when pL ^> 1 the zero mode effectively disappears. Thus 
standard gravity only appears at energies below the compositeness scale L -1 and is 
seen to emerge from the CFT, whereas above the compositeness scale (or r ^ L) 
the dominant CFT contribution changes Newton's law to a different power law (J84|) . 
which can be fractional and even irrational. 

The cancellation of the leading 1/p 2 term can also be seen in coordinate space for 
special values of a where the sum over Kaluza-Klein modes can be done analytically. 
Consider the case of z/_ = —1/2 or a = 3/8 (recall that a > 1/4 corresponds to 
IR-brane localization). The Green's function is: 



G p (z , z ) = -- cothQ^i) = — - 



P 



2 £ 



1 



^ p 2 + m 2 



(90) 



where m n = nn/zx (note that this, in fact, is the Green's function for an even 
field on a flat S 1 /Z 2 orbifold of size Z\, evaluated at one of the fixed points). It is 
straightforward to check that the 1/p 2 pole is cancelled for pz\ ^> 1. Alternatively 
the gravitational potential can be calculated exactly using the expressions ()D.20|) and 
()D.23j) in Appendix D. One obtains: 



V(r) 
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(91) 



Note that the Fourier transform of V is proportional to the Green's function ()90|) . 
For distances r < z\ ~ L the series expansion gives: 



V(r) 



fi 2zi 
' M 2 p Tir 2 



+ ... 



(92) 



This agrees with the gravitational potential (J84|) for r < L = 1/ (kAi). 

For a < 1/4, one can also analytically verify that the 1/r term does not cancel 
when z/_ = 1/2 (or a = 1/8). This is obvious because in this case the exact Green's 
function is: ^ ^ 

Gp{z , z ) = — coth(^i) , (93) 

and the extra —1/p 2 term compared to (J9*U|) . corresponding to the 1/r term in the 
potential, is no longer cancelled at high energies. 

These two simple examples enlighten our gravitational potential calculation for 
a > 1/4, and reveal that when the graviton is localized on the IR brane the usual 
Newton's law of gravity only emerges as a low-energy phenomenon. Remarkably, as 
we will see in the next section, these properties can also be described purely in a 4D 
dual interpretation. 
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5 The holographic interpretation 



Motivated by the string theory AdS/CFT correspondence, bulk theories in a slice 
of AdSs can be given a holographic interpretation as dual to a CFT (at large- N 
and large 't Hooft coupling) with conformal invariance broken in the IR, coupled to 
gravity and possibly other fields (SHEIIZI- As opposed to the string-theory AdS/CFT, 
the UV boundary value of bulk fields are not only sources of operators in the dual 
CFT but acquire their own dynamics due to the presence of a UV cutoff. 

The nature of the holographic interpretation of bulk theories in a slice of AdSs 
is necessarily speculative; however, it passes some quantitative tests, despite the fact 
that one usually does not know what the dual CFT is, or even whether it actually 
exists! Nevertheless, in most cases one is able to give a reasonably convincing picture 
of the dual 4D strong dynamics, which may be useful as a guide to finding a theory 
with the desired features. With these caveats in mind, in this section we attempt 
a quantitative dual 4D description of the mass-deformed gravity in a slice of AdS5 
studied in the previous sections. 

The study of the 4D dual begins by recalling that the UV-boundary values of 
the bulk fields are sources of operators in the dual CFT. The classical bulk action, 
evaluated on solutions of the bulk equations of motion with arbitrary values on the UV 
brane (and obeying the proper boundary conditions on the IR brane) is interpreted 
by the stringy AdS/CFT correspondence as being equal to the generating functional 
of connected Green's functions of the 4D dual CFT. To continue, we begin with the 
solution of the bulk equation which obeys the boundary condition on the IR brane: 



H(p, 



H(p)A- 2 (z) ( J u±Tl (iq) -Y u±Tl (iq) 



K± (iqi) 



(94) 



where v± = ±(4a± — 1), q 
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—m 2 is the mass- 



shell condition, and H(p,z) is the 4D Fourier transform of h(x,z). Note that we 
will omit tensor indices, as we are only concerned here with the transverse-traceless 
components of the metric perturbation, which obey a scalar equation. The bulk 
action, which generates Green's functions in the dual is: 

M 3 f d 4 p 



S, 



bulk 



A 3 H{p, z)(H'(-p, z) - 4aAkH(-p, z)) 



4 J (2tt) 4 

evaluated on the bulk solution (|94|). This can be rewritten as: 

M 3 k f d 4 p 
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s, 



balk 



(2vr)< 



F(q ,qi) H(p)H(-p) 



(95) 



(96) 



where: 



=R9o 



Ji 



TlW0. 



- Y, 



Y u± (iqi) 



Ju ± (iqo) ~ Y u± (iq 



Ju± (iqi 



(97) 
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The dual theory two point function of the operator O sourced by the bulk field 
H, (00)(p) = f d 4 x e~ ip ' x (TO(x)O(0)) , is contained — up to local counterterms 
which we will discuss later — in the second derivative, £(p), of Sbuik with respect to 
the boundary value of the metric perturbation A^h. The correlator is, in various 
equivalent forms to be used later: 

E(p) = / d 4 xe~ ip - x ^ 52Sbulk ^ 

J 5(A 2 h(x,z ))5(A 2 h(0,zo)) 

M\ 3 k 4 . . Jv±(iqo)Yv±(iqi) -Y v± (iqo)J v ±(iqi) 
-r TT ! : m 



k) 2 J v± ^i(iqo)Y v± (iqi) -Y Vi ^\(iqo)J v± (iqi) 

M\ 3 fc 4 g (/,(go)/C(gi) - I v (qi)K v (q )) 



k J 2 I UTl (q )K u (q 1 ) + I v {q\)K u ^ 1 {q Q ) 

M3 1 (98) 
2Aq G p (zq, zq) 

where G p is the boundary-to-boundary propagator of (jHBj) divided by A^ and with k 
replaced by A k. 

The behavior of S(p) can be studied for momenta p such that kAo ^> p ^> 
kAi, or equivalently, ^ ^ > 1. In this energy regime, the effects of the 

conformal symmetry breaking (i.e., the IR brane) are completely negligible. The 
leading nonanalytic piece in S(p) is then interpreted, by "matching" to the string 
AdS/CFT correspondence in the A — > oo limit, as due to the strong dynamics of 
the dual CFT above the scale of conformal symmetry breaking. On the other hand, 
the analytic pieces in the correlator, which, in string AdS/CFT, are subtracted away 
by adding appropriate counterterms, are now interpreted as kinetic (and higher- 
derivative terms) of the dynamical source field in the holographic dual. 



5.1 a_ branch holography 

We will first consider the ct_ branch of our solution, which is the one continuously 
connected to the a = RS value (see Figure |TJ). For | > ct_ > 0, or — 1 < z/_ < 1, 
we find: 8 



M\ 3 k 4 



r(-u) 



(99) 



while for all other values on the a_-branch (a_ < or v = z/_ > 1): 

3 

2 



E(p) 



M 
T 



k 4 \q% 



2u 



2u] 



l + ... + cq l J +q z 



J2u+2 



-V) 



vV(u) 



+ 



(100) 



8 Notice that the expression (|99J1 does not contradict l|89Jl since when < £ < 1 (or — 1 < v- < 0) 
the nonanalytic term dominates in (|89|l . so that the constant term disappears in the inverse. Also, 
for 2/_ = —1/2 there is no analytic term in the series expansion in this regime. 
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Figure 2: The scaling dimension Aq plotted as a function of the boundary mass 
parameter a. The source coupling to the CFT is irrelevant for a < 1/4 and a > 3/4, 
marginal for a = 1/4 and a = 3/4 and relevant for 1/4 < a < 3/4. Note also that 
the RS value at a = 0, and the specular RS value at a = 1 both have A© = 4. 



where [2i/] denotes the largest integer smaller than 2v. In each case (199(1 . ()100J) . 
we have included the leading analytic piece as well as all terms up to the leading 
nonanalytic piece, q^ v+2 - The power of go i n the nonanalytic piece indicates that the 
scaling dimension Ae> of the operator O — the energy momentum tensor of the dual 
theory — sourced by the metric perturbation h is: 

A = 3 + z/ = 4-4a_ , (101) 

on the a_ branch. The leading analytic piece in (|99j) . (|1UU|) indicates that there is a 
kinetic term for the metric perturbation in the dual theory. 

Thus, the holographic description of this branch is that of a metric fluctuation 
h^ u coupled to T^ V FT of scaling dimension 4 — 4a„ . Note the unusual fact that the 
energy momentum tensor of the CFT has an anomalous dimension. This, however, 
is required if the 4D dual is to evade the Weinberg- Witten theorem — which assumes 
Poincare invariance, broken here by the presence of a nontrivial background metric, 
as in theories of induced gravity. Notice also, from (|101|) . that the scaling dimension 
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of T£ FT can be as low as 2, for 1/2 > a_ > 0, as can be seen in Fig. El A scaling 
behavior of the CFT with such a large anomalous dimension should persist no matter 
how small the breaking of Poincare invariance (or, equivalently, the deviation of the 
metric background from Minkowskian) ! Leaving aside the issue of existence CFTs 
with such behavior, we continue with our attempt at giving a (semi-) quantitative 
picture of the dual dynamics. 

5.1.1 The dual theory and its dynamics 

The Lagrangian of our dual theory is, then, at a UV scale ~ k, with a canonically 
normalized metric perturbation: 

£uv = e\ K P Uh^ + ^ h, p T^ FT + ^f- h^TZuer + £cft , (102) 

where e u = sign v and Xuv — |^| 1//2 (M/fc) _3//2 . We have included the coupling to 
observable matter fields (UV-brane localized in the gravity dual). 9 

From eqn. (jl(J2|) . taking into account the anomalous scaling dimension of Tqft 
from (jlOlj) . we conclude that the coupling of the metric perturbation to the CFT en- 
ergy momentum tensor is relevant for a_ > 1/4, marginal if a_ = 1/4, and irrelevant 
for a_ < 1/4. Introducing a renormalization scale //, the dimensionless coupling is 
then A(/i) = (^/k) 1 ~ 4a " Xjjy and satisfies the RGE: 

^ = -(4a_-l)A + ... (103) 
dfx 

where the first term is a result of dimensional analysis and higher order terms due to 
the CFT's interactions have been neglected. 

If A is relevant, i.e. a_ > 1/4, the solution of the RGE: 

(, \ 4a--l 
, (104) 
m IR J 

indicates that the coupling of h^ v to the CFT is enhanced in the IR. The conformal 
invariance is broken at the IR scale tuir and we expect that integrating out the CFT 
dynamics at the IR scale will induce a kinetic term for h^. From the CFT point of 
view the dynamics is both strong and unknown; however, it is clear (see, e.g., [3*U] ) 
that producing a kinetic term requires two insertions of the dimensionless coupling 
Xir, as indicated in eqn. (J106J) below. 

Clearly, in the weakly coupled gravity dual we can directly compute this contri- 
bution by calculating the two point function E in the IR limit p <C kA\. The pure 

9 Note that in the case of irrelevant coupling, the sign of the kinetic term for h^v is the proper 
one, as t v — 1 for a_ < 1/4; while it has the wrong sign in the case of relevant coupling; we will see 
below that the leading contribution to the h^ p kinetic term of the right sign arises from IR physics 
not accounted for in (|102|l . 
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IR contribution is obtained by subtracting the analytic piece of £ that arose in the 
limit p ^> kA±, see (jSHJ)- Thus, expanding (J9*%|) for small q and q\ and subtracting 
the analytic term already accounted for in (|99|) (and in the kinetic term in (|102p ) 
leads to the pure IR contribution to the correlator: 

where A\ = mm/k. Thus, the effective Lagrangian describing the long-wavelength 
fluctuations of h^ u and its coupling to the observable matter sector is given by: 



IR 



MYXj R 



k j v 



- A K p Uh^ + ^Lh^TZtter ■ (106) 



In the relevant case, where —l<v_< 0, our gravity dual calculation shows that the 
IR contribution to the kinetic term for h w has the correct sign and dominates over 
the ghost-like e u contribution. Hence, using (jl(J4j) and canonically normalising (|1U6J1 . 
the coupling of the observable matter to gravity at scales below mm is: 



k\u-\ Xuv V k\u-\ \m IR J 
Upon identifying the ratio of UV to IR scales with the warp factor, 

— = e nkR , (108) 



m IR 



Eq. p()7j) gives Mp = e 27rfc_R(4«_ i)^ j n a g reemen t with the gravity dual result 

dHDJ)- 

Consider next the irrelevant case with u_ > 0. In this case the coupling of 
the metric perturbation to the CFT is irrelevant so that the induced contribution 
proportional to \m in (|106J) is negligible. Thus, in ()1U6|) the leading term proportional 
to e v dominates (and has the right sign!), leading to M P = k/Xuv, which precisely 
equals the last line in Eq. (|50|) as well as the RS result. 

Finally, consider the case where the coupling of CFT to the background metric 
is marginal (z/_ = 0). We have to take into account higher order terms in ([103)1 . 
which we can calculate using the weakly coupled gravity description. In the IR limit 
p kA\ we obtain: 

f M\ 3 k 4 ( 2 A 



S(p)« = - ( T j ~ U log + ... ) , ( 1 09) 



so that the IR Lagrangian becomes: 

3 / a \ 2 



(M\ (A \ 1 1 
TJ l h ^ UhW + lK P TZtter , (HO) 
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where Aq/A\ = k/rrij^. Thus canonically normalising the kinetic term and using 
(jl()8j) leads to a Planck mass: 



This again agrees with the corresponding result in (|80l) . 

Thus, we have a rather unusual "theory," particularly in the case of a relevant 
coupling of the CFT to gravity. We have two important hierarchical scales of na- 
ture, for defmiteness take them to be meV (m IR ) and TeV (k), as in our discussion 
after Eq. (|SUJ). The infrared scale is presumably determined by some dynamical 
mechanism — unspecified both in the 5D gravity and the 4D CFT descriptions — from 
the UV scale. The TeV scale is the cutoff of the theory, where a more fundamen- 
tal description takes over. Naturally, in the UV theory, a Newton constant Gn of 
order TeV -2 is expected (the first term in (jl(J2|) . possibly including additional bare 
contributions). Along the RG flow to lower and lower scales, the hidden CFT sector 
becomes stronger and remains so until, at the meV scale, conformal invariance is 
broken in the strongly coupled CFT. Thus, the meV scale "broken CFT" induces a 
Newton constant. It is strong enough that despite the fact that it operates at meV 
scales, the induced Planck scale is hierarchically larger: Mp 3> TeV 2 ! In other words, 
gravity is so weak because of the strength of the hidden CFT over a large interval of 
scales. 

While this "scenario" sounds really unusual, and we are not aware of a CFT with 
the desired properties, it is not so difficult to come up with a weak coupling — in fact, 
free field theory — model of this phenomenon. Consider N free fields (of whatever 
nature, so long as they induce the correct sign Newton constant), coupled only to 
gravity, with a characteristic mass scale ~meV. Above the meV scale, this hidden 
theory is conformal, which protects the UV modes of these fields from generating 
a Gn. But conformal invariance is broken at the meV scale and so one expects a 
contribution to Mp which will be of order iVxmeV 2 . Now if we take N = Mp/meV 2 , 
then this clearly dominates the TeV 2 contribution of ordinary massive matter. 

As far as the evolution of the universe, our picture would predict that the strength 
of gravity should change from TeV to Mp during cosmological evolution; note that 
this does not affect BBN as the transition to "normal" strength gravity occurs before 
nucleosynthesis, when the Hubble size is of order meV^ 1 . There may, however, be 
relevant consequences for the physics of the earlier universe, as for instance inflation 
and baryogenesis. 10 

5.1.2 The gravitational potential in the dual theory 

Let us now describe the leading correction to Newton's law at intermediate distances, 
r < l/(kAi), from the point of view of the dual interpretation. In the irrelevant case 

10 It should be clear that we have nothing to say about the cosmological constant. 




(Ill) 
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Figure 3: The Feynman diagrams in the 4D dual theory responsible for the gravita- 
tional potential corrections (which need to be summed up for the case of a relevant 
coupling). The source field h^ u , indicated by a wavy line, interacts with the CFT 
contribution, indicated by the blob. 



(u- > 0), the coupling of matter to the CFT can be treated perturbatively and 
the leading correction arises from a single insertion of the CFT correlator and two 
insertions of the source field, as in the RS case (see Figure EH). The leading and first 
subleading contribution to the Newton potential is: 

>?uv fd 3 P Jpx (I \ 2 uv (00){p)\ 



V(r) = ^ e 



k 2 J 2n 2 \p 2 k 2 p 4 J 

- n X uv f d3p r ip - x (rr 2 n 2u ~ 2 k- 2u ^ (^^?) 

* k 2 J 2ix 2 6 \ P P k 2*T(i/) J ' { } 

where: 

is the nonanalytic piece of S(p) of ()99jh interpreted as the CFT correlator at the 
relevant energy scale, and Xuv is given after (jl02|) . A Fourier transform is performed 
by using the properly regulated and normalized 3d Fourier transform of p a , which 
is jSU: 

P e^p«=r-) - . (114) 



2tt 2 r \rj 2V5Fr(-f) 

After using ()8U|h and various gamma-function identities, 11 we find precisely our result 
from the gravity calculation (|%3*j). 

Consider next the correction, at r < l/(kA\), for the case when the interaction 
with the CFT is relevant (z/_ < 0). Then, we have to sum the chain of bubble graphs 
as indicated below (recall e v = — 1 now): 



V{r) = -fi 



X 2 UV f d 3 p e ip - x 



k 2 J 2n 2 p 2 



?vv (QQ)(P) + 6v Kv ( (QQ)(P) \ 
k 2 p 2 v k 4 \ p 2 J 



x h_ f 1 „ ffl ..„,_ 



II uv / £_ P *P- X ~ — /; / — p z P' x fll^l 

^ k 2 ] 2^ 6 p2 _ %v {0 0)(p) ~ (00>(p)' ' ' 



11 A useful identity is: T(2x) — 2 T(x)T(x + |). Note also that a factor of l/(47r) has been 
absorbed in the definition of V. 
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where we notice that for the distance scales of interest the CFT correlator dominates 
over p 2 in the denominator, as appropriate for a relevant coupling. Finally, computing 
the Fourier transform as before, and using (jSUj) we again recover precisely the leading 
term of the potential from Eq. (|84p from the gravity side. 



5.2 a + branch holography 

Let us now consider v — v + — 4a + — 1 > 1. In this case the graviton is always 
localized on the IR brane. We find, for A k ^> p Aik, taking the upper sign in 



M \ 3 

~ ( T J fc 4 (116) 



x 



I i\ , 2 1 , [2i/]-2 . 2^-2 T( 2 ~ ^) 



4(z/-2) " u ' " u 2 2 - 2 r(z/-l)_ 

indicating that now that the scaling dimension of the operator O is 

A = v+ 1 = 4a+ . (117) 

Thus when u + > 2 the source coupling to the CFT is irrelevant, marginal for v + = 2, 
while for 1 < u + < 2 the coupling is relevant. This behavior is plotted in Fig. EJ 

At low energies pL < 1, on the other hand, the series expansion of the correlator 
is given by: 12 



Z(p) IR ~- — k 



.\r 3 



A 



118) 



To obtain (|118j) . we first took the large- A limit in (J98|) . obtaining a power series in q 
(of which we kept only the leading terms above), with (^-dependent coefficients, and 
subsequently expanded these coefficients for small q±; a more formal way of obtaining 
this is by multiplying the entire correlator by Aq U ~ 2 and taking the large Aq limit (see 
footnote). Remarkably we see that the correlator now has a pole at p 2 = 0. Thus at 
low energies we can interpret the massless graviton to be predominantly a composite 
of the CFT. This is in contrast to the a_ branch where no such pole exists. Similar 
pole structures have also been identified in Refs. [HI HH E2] ■ 

12 Note that (|1 18|) does not contradict (j88(l . which asserts, instead, that ~ G(p) ~ 1/p 2 

for p < L _1 . In l|118|l the pole in E at p 2 — appears only if the UV cutoff is taken larger than 
the AdS curvature scale; instead in (|88|l . finite cutoff effects (i.e. mixing with the source) wash 
away the pole. Note that if we formally took Ao — * oo (the UV cutoff of the dual theory), after 
multiplication by A 2 ^ -2 to isolate the leading nonanalytic piece in (|If 6f) . the persistence of the pole 
in (|II8|) would run afoul of the Weinberg- Witten theorem. Our dual field theory i|119|) . however, is 
coupled to gravity perturbed by the graviton mass, hence the cutoff cannot be taken large (and the 
source decoupled), due to the strong coupling problem of massive gravity |26) . 
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In addition, we immediately see from ()118)) that the leading analytic piece is 
a constant, and corresponds to a mass term of order the curvature scale for the 
source field h^ v . The interpretation of this fact in the dual 4D theory is that the 
CFT generates a mass for the source so that it decouples at low energy and the 
propagating mode is predominantly the composite graviton (see also the calculation 
of the long-distance gravitational potential (|121|l ). 

The analytic terms of (|118|) can be used to obtain the long-distance Lagrangian: 

C IR = \ h w {U - m 2 h )h^ + | h w T£ p FT + | h w TZtter + ^cft , (119) 

where x — iy+ ~ 2) 1 ^ 2 (M/k)^ 3 ^ 2 and m\ = 4(z/ — l){y — 2)k 2 . If we now write the 
small-momentum expansion of the correlator as: 

(OO) ~ (Mkf 16i/+(i/ + - l) 2 ^f+4 , (120) 

where Aq — 1, then the leading contribution to the gravitational potential at large 
distances is given by: 

V(r) - -u£ f^e'^™ 
V[r) ~ ^k 2 J 2n 2 k 2 m\ ' 



X 4 M 3 ___ ^ 2a2u+ fd 3 p oip . x l 



m\ k J 2ir z p A 



M 2 r 



(121) 



where the source propagator has been approximated by l/(p 2 + mf l ) — l/ m /i and the 
Planck mass is given by (using the values for and x given after (|119J) ): 

This agrees with the Planck mass formula ()8()|) for a + > 1/2. Note also that further 
insertions of (OO) in ()121|) are negligible at large distances. 

In fact, the result for the leading long-distance contribution to the potential (|121|) 
indicates that at r > L, we can describe long-range physics as due to the exchange 
of a massless spin-2 field, — mostly a composite of the CFT — coupled directly to 
matter, thus forgoing the discussion of the massive source in (|119|) . In this case the 
Lagrangian is given by 

1 / M\ 3 k 2 

ClR * ~ A (j) - K 2V K p UW» + K p TZtter , (123) 
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where canonically normalising the kinetic term leads to the Planck mass (jSOjh The IR 
lagrangian (j!23j) can, equivalently, be obtained by directly considering the IR limit 
(qo -C 1, qi <C 1) of the correlator of 





v 3 r 


(f 


) ¥ 



(124) 



and interpreting the leading A[ 2v+ term as the kinetic term of the interpolating 
long-distance field h^ p in (jl23j) . 

Continuing on to intermediate energies L™ 1 <C p <C k, we see immediately from 
(jll6|) that there is no longer any pole, since we are now above the compositeness 
scale. This is also consistent with the fact that for pL ^> 1 there is no longer any 
1/p 2 term in the Green's function (|89j) . The analytic terms are identical to those at 
low energies, so that the source remains massive. Hence, in this regime the dominant 
contribution to the gravitational potential arises from the CFT, which corresponds 
to the nonanalytic term in f)116j) . 

When v + > 2 the source coupling to the CFT is irrelevant and the gravitational 
potential follows from the coupling to the CFT as depicted in FigEl From the analytic 
terms of we obtain the UV Lagrangian: 

Cuv = -g hfj, p (D - m 2 h )h^ p + j- h w T^f FT + j- hp, p T^ atter + Ccft , (125) 

where x = { v + — 2) 1 / 2 (M/A;)~ 3 / 2 and w? h = A{y — \){y — 2)k 2 . Since there is no longer 
any massless pole, the leading contribution to the potential is given by: 

V(r) ~ f d *P c iv* X 2 (00)(P) 

{ } ~ ^k 2 J 2vr 2 k 2 m\ 1 

X 4 (M\ 3 f d*p / p y»-z AT{2-u) 

= ^<\JJ J^~ 2 ^) r(i/-i) ' (126) 

where (OO) is the nonanalytic part of given by: 

Performing the Fourier transform leads precisely to the result derived purely on the 
gravity side (|84|. 

When 1 < v + < 2 the source coupling to the CFT is relevant but the nonanalytic 
term is still subdominant compared to the leading mass term in f)116j) . In this case no 
summation is needed beyond the leading CFT correction depicted in Fig. Eland so the 
contribution to the potential is identical to that obtained in (j!26|) . The corresponding 
Fourier transform then leads to the same expression (|84j) . 
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6 Discussion and conclusion 



We have seen that the graviton zero mode can be smoothly deformed away from 
the Planck brane. This deformation requires modifying the bulk covariant theory at 
quadratic order by introducing bulk and boundary mass terms. However, a massless 
mode only occurs for a special choice of the bulk and boundary masses. This is an 
additional tuning beyond the usual tuning of bulk and brane cosmological constants 
in the RS model. As in the RS case this tuning may be realized as the result of bulk 
supersymmetry EH] - At the linearized level there is a 4D general covariance, which 
is consistent with the fact that there is a massless tensor mode. In addition there is 
also a massless vector mode with a 4D U(l) gauge symmetry. 

However, general relativity is an inherently nonlinear theory, and it is apparent 
that higher order nonlinear interactions will spoil this symmetry. Without further 
modification the gravity in our model is different from the full nonlinear Einstein 
theory. This situation may be remedied by modifying our model at the nonlinear 
level with the introduction of nonlinear terms in the bulk and brane, in order to at 
least preserve the 4D general covariance. Thus we expect our zero mode to remain 
massless in the nonlinear theory, although this analysis is beyond the scope of the 
present paper. 

Nonetheless it is already interesting that a smooth deformation exists at quadratic 
order without the presence of ghosts. In particular the scalar sector is trivially zero 
because this is the only solution consistent with the bulk and boundary equations. 
Clearly this is due to the fact that we are working only to quadratic order, and the 
scalar modes can possibly appear at the nonlinear level. Scalar modes may also arise 
when matter is added on the brane. On the phenomenological side, they are certainly 
needed to reproduce the correct gravitational law if the stress-energy tensor of the 
matter fields is not traceless. A similar situation takes place in the usual RS case. In 
the compact version with a stabilized radion there are no massless scalar excitations. 
However, a massless scalar mode (most easily interpreted as the brane bending mode) 
arises when matter is present on the brane, and allows for the recovery of standard 
4d gravity at large scales j^lj. A similar analysis should also be carried out in the 
set-up we have discussed here. 

By the AdS/CFT correpondence there is an interesting 4D dual interpretation 
of our model, especially in the case when the graviton zero mode is localized on the 
IR brane. This is because zero modes localized on the IR brane correspond to CFT 
bound states and therefore the dual CFT interpretation would correspond to gravity 
emerging from the strongly coupled gauge theory. In this model of emergent gravity 
the UV theory is a gauge (string) theory at the TeV scale, and the graviton is a 
composite particle which can be associated with the millimeter scale. Thus gravity 
emerges as a low energy phenomenon in the IR. This is different from the conventional 
viewpoint that gravity is a fundamental degree of freedom in the UV theory, and our 
model is the first step in constructing and understanding this novel possibility. 



33 



Acknowledgments 



We would like to thank Joel Giedt and Alex Pomarol for useful discussions. The work 
of T.G. and M.P. was supported in part by a Department of Energy grant DE-FG02- 
94ER40823 at the University of Minnesota. T.G. is also supported by a grant from 
the Office of the Dean of the Graduate School of the University of Minnesota, and 
an award from Research Corporation. E.P. acknowledges the support of the National 
Science and Engineering Research Council of Canada. T.G. also acknowledges the 
Aspen Center for Physics where part of this work was completed. 



A Bulk gravity action to quadratic order 

We will present the expansion of the bulk action (fTtjj) around the background RS 
solution, to^ quadratic order in the perturbation h MN , where gMN = A 2 r] MN + 1imn = 
A 2 (r]MN + hMN), and M, N — 0, ... ,3,5. The Lagrangian density to quadratic order 
in h M N is given by: 



C 5 [A 2 V + h; M] 
M 3 A 3 



h MN d K d K h MN - 2h MN d N d K h MK + h MN d M d N h + ~d K h MN d K h MN 



+ d N hd M h 



MN _ \d K hd K h - U K h MN d M h KN - d M h MN d K h KN 

+ H d M d N h MN - \h d K d K h + + 2^ + A 2 ^j (h MN h MN - h 2 

i ~~ \ A' 2 (~ ~ 
-hh'\ +12— ( h A5 h A5 -^hh 



A 



+ 4^ I h MN U 



MN 



A 2 



1- 
— / 

2 



,55 



- A 2 k 2 (ah M N~h MN + b~h 2 ) 



(A.l) 



where h = h^, prime (') denotes d$, and indices are raised and lowered with t/mn- The 
last term in (jA.l|) is the contribution from the bulk mass term. The corresponding 
equation of motion arising from £ 5 is given by: 



M 3 A S 



- d (h MN - rj MN h) - r] MN d A d B h + d M d h AN + d N d h AM - d M d N h 



A' f~ ~ ~ ~ ~ \ A' 2 ~ 

— 3— ( h MN — rjMNh + 2r]MNd A h — dj^hm — <9tv^a/5 j — 12-^-^mtv^55 

-2A 5 (A + 6k 2 M 3 )(h M N - -VM N h) + 4A 5 k 2 M 3 (ah MN + brjMN h) = , (A.2) 

where d 2 = □ + d 2 . Note that the R-S solution jT] requires that A = — 6k 2 M 3 , and 
the term involving the cosmological constant in the last line of (|A.2|) vanishes. If we 
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consider only the tensor fluctuations Hmn = as defined in (jHJ) then we recover the 
equation of motion ()18j) . The remaining equations of motion for the vector (|23|) . (|24jl 
and scalar modes ([30 )1 - ([33 )1 follow from the /i5 and 55 components of (|A.2|) . A similar 
expansion to quadratic order in the metric perturbation, but without the bulk mass 
term, has also been performed in Ref. [2*Tj . 



B Quadratic action of the graviton zero mode 

We have seen that by appropriately tuning the bulk and brane mass parameters, there 
is a zero mode both in the tensor and vector sectors. The presence of these zero modes 
signal some symmetries of the starting action. In particular, the massless tensor mode 
signals a 4D general coordinate invariance, while the massless vector mode signals a 
4D gauge invariance. Clearly, these symmetries only occur at the linearized level in 
the perturbations. A study beyond this order would require considering the inclusion 
of higher order terms in the original action, and it is beyond our current aims. 

It is instructive to compute the quadratic action for the tensor zero mode, and 
to explicitly show how the mass term cancels. This calculation also reveals how to 
canonically normalize the graviton in the 4D theory, and therefore how to properly 
define the four dimensional Planck mass. Assuming the Fierz-Pauli choice b = —a, 
and Pi = —aii, with —ol\ = a = a for the bulk/brane mass terms, respectively, the 
total action up to second order in the tensor perturbations gives: 



(2) 



tensor 



d 5 xA 3 



1- - A' ~ ~ / A' 2 
— h'^ v h' h' - 3 — + 3 k 2 A 2 + a k 2 A 2 ) hTh 

4 A * \ A 



± / <Fx 



A ! k I ^-a)/^V 



(B.3) 



where +/— refers to the UV/IR brane, respectively, and the h^ u spacetime indices 
are raised with the (inverse) Minkowski metric r]^ . The first term, 



a/-?4 Ra 



(B.4) 



has the tensorial structure of the quadratic 4D Einstein-Hilbert term for a transverse- 
traceless h^, except that h^ v still depends on the fifth coordinate z. One can verify 
that the action (jB.3|) reproduces the tensor mode bulk and brane equations, (JTHj) 
and (|3Hjl. respectively. 

To proceed further, we decompose h^ u into the eigenmodes (fE?|) and ()2(J|l. and 
integrate over the compact coordinate. Hermiticity of the action ensures that eigen- 
modes with different mass eigenvalues are decoupled, so that one is left with an 
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infinite sum over the decoupled actions S^^, for each four dimensional mode H$ . 
The last two lines of (|B.3|) combine to form the mass term for any given mode. In 
particular, let us consider the zero mode: 



h 



C x A{z)-^-^ Hf){x) 



(B.5) 



as given by (|19j) . where we only include the part which is continuously connected to 
the RS graviton. In this case, the mass terms of (jB.3|) combine to give: 



C^kM 3 
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d 4 x H^H^J. (B.6) 



Indeed we see that the bulk and brane contributions cancel when (I64|) is imposed for 
a = a_, resulting in a massless tensor perturbation. Similarly, if we choose the C2 
part of (pPIJ) then we obtain a massless tensor perturbation for a = a+. 
The quadratic 4D action for the zero mode then simply becomes: 



c(2)(0) 
scalar 



Zl 



,2(l-4a) 



A(z 



,2(l-4a) 



] J d 4 x^f,R, , (B.7) 



2A;(4a 

where , i?4 now refer to the standard 4D metric = rj^ + H^){x). By adding 
a source term on the brane, it is straightforward to see that the coupling of H$) to 
brane fields is set by the 4D Planck mass J7SJ). 



C Quadratic action of the scalar modes 

The equations for the scalar perturbations can be directly obtained from the second 
order action in the scalar perturbation. The derivation of this action is quite involved 
(we extend the computation of Ref. [33], performed for the covariant case), but the 
result is relatively simple. The computation of the action supports the choice of the 
(generalized) Fierz-Pauli mass term for the perturbations. Indeed, the expansion of 
the mass terms gives the following higher derivative kinetic terms: 

^cL ^4 A; 2 J d 5 xA 5 (z)(a + b)(nE) 2 + J2^ k J <?x A\z t )(a, + A) (□£)■ , 

(C.8) 

where the sum over % refers to the two boundary branes. As in the 4D case, the Fierz- 
Pauli choice (a+b = ai+f3i = 0) eliminates these pathological higher derivative terms 
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from the action [26J. Employing this choice, and relating the brane mass coefficients 
as «o,i — , the action for the scalar perturbations is: 



(2) 

scalar 



d 5 xA s 



X 



6$D + t) + 12^ 2 + 24^' (2 $ - J) + 12^- ( J 2 + 8 
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± / d 4 xA 3 <j24^ 2 + 6^nc + 3^cnc 



+ akA 



4 8(^ + ^c) +24^ + ^C)^ 



(C.9) 



where +/— refers to the UV/IR brane, respectively. Note that only gauge invariant 
combinations xp, <fr, and£ appear when there are no bulk or boundary masses. This 
is due to the fact that the usual action for gravity is general coordinate invariant. 
This symmetry has been made manifest in (|C9|) . by rewriting a total derivative in 
the bulk as a boundary term, which then produces a sum of two brane terms. For 
instance, this is the origin of an E" term which is not present in the original action, 
but which is needed to produce the ip 2 term in the bulk. Clearly, this procedure does 
not change the action, but allows one to write it in a more compact and manifestly 
covariant form. The terms proportional to the bulk/brane masses are instead not 
general coordinate invariant, and for this reason they cannot be rewritten in terms 
of gauge invariant quantities only. 

To obtain the equations of motion from (jC.9|) . it is convenient to work with the 
original total derivative bulk terms, and compute the Euler-Lagrange equations by 
the usual procedure. In fact the simplest way to proceed is to rewrite (jC.9|) as a bulk 
action by promoting ( to a bulk field, which evaluates to Q at the two boundaries, 
and then vary this 5D action. In general when we vary terms containing /' (where 
/ denotes any of the scalar perturbations), we produce terms proportional to Sf. 
These terms are dealt with by integrating by parts the variation of the action: 



d 5 x5f [...]= / d 5 x{5f [...]}'- / d 5 x5f [...] 



(C.IO) 



The last term in (|C.10|) enters in the usual Euler-Lagrange equations, while the total 
derivative is usually assumed to vanish when evaluated at infinity. However, in the 
presence of branes, this term must also separately vanish on shell, for each brane. It 
is precisely this requirement that leads to the boundary conditions for the bulk fields. 
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Using this procedure, the bulk/brane equations of motion can be readily com- 
puted. There are four bulk equations: 
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(C.ll) 



which are obtained by varying B ,E , , and ip , respectively (since C is defined only 
on the two boundaries). 

For the boundary conditions one would naively expect five equations. However, 
one can verify that once the boundary term in (|C.9J) is rewritten as a bulk term, <$' 
and B' do not appear in the total action. Hence, there are no boundary conditions 
arising from the variation of the action with respect to and B. This leaves only 
three boundary conditions which follow from varying E ,ib , and ( in the action (jC.9|) . 
These are respectively given by the following equations evaluated on the two branes: 



□ 



A' ( A' 

<//-— c(>-AakAU + — C 







(C.12) 



A' f A' \ 1 

ip' - —cf>-4:akAU + —() + -D[E' - B - ( - 4akAE] = , (C.13) 



A , A' ( , A' \ □ 



. 



(C.14) 



For the massless modes (i.e, imposing □ = 0) the solution of the above equations 
is , = = C = O. The remaining combination B — E' is undetermined and one 
can verify that when the other modes are absent, this combination gives a vanishing 
contribution to the action (jC.9|) . In the covariant case, we have seen that out of the 
original scalar modes ib, 0, E, B, £ only ip, 0, ( appear in the quadratic action for 
the perturbations. This is a consequence of 5D general coordinate invariance, which 
guarantees that two modes are not present in the action at all orders. In the present 
case, higher order terms could make the remaining modes dynamical. However, to 
study in a consistent way the non-covariant theory beyond quadratic order in the ac- 
tion requires introducing higher-order terms in addition to the quadratic bulk/brane 
mass terms and (|5U|) . This is beyond the scope of the present analysis. 
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For the massive modes, one can check that the bulk equations (jC.llJ) are equiv- 
alent to the equations ((HP)) - (jHSJ)- Also, the boundary equations (|C.12J) and (|C.13J) 
are equivalent (for the massive modes) to the two junction conditions (|60|). and (|61|) . 
However Eq. (|C.14|) is a third independent boundary condition, which as discussed 
in the main text, cannot be obtained through the usual junction/Israel condition 
procedure. Eliminating ( in (jC.14j) through the other two boundary conditions gives: 

Using the third equation of (jC.llJ) . this equation simplifies to: 13 

a(4iP + DE) = . (C.16) 

This boundary condition vanishes identically when a = 0. This shows that the "ad- 
ditional" boundary condition (|C.14J) is redundant in the standard case without bulk 
or boundary mass terms. Instead when a ^ we obtain an independent boundary 
condition, which appears as Eq. in Section 3.1. 



D Derivation of the gravitational potential 

We will now derive the leading terms in the expressions (|83j) - (j84j) . We start from the 
action for the Kaluza-Klein modes of the 5d tensor mode, coupled to a conserved 
matter source on the UV brane. The action can be written in the form: 

S = c n I d 4 x HV»" (□ - m 2 n ) H$ + ±Z n (0) f d 4 x H$T»" , (D.17) 

where denotes the energy-momentum tensor of the matter source, and the vari- 
ables of the tensor mode wavefunctions are separated as: 

h flv (x,z)=Y,H$(x)Z n (z) . (D.18) 

n 

The coefficients c n are given by: 

m 3 r zi 

c n = —J o dzA 3 Zl(z) . (D.19) 

Correspondingly, the gravitational potential generated by a static mass \i on the 
UV brane, measured at the distance r from \x , is given by: 



13 Note that there is no problem with combining bulk and boundary equations. Indeed, the value 
of any function f(zi) in the boundary equations is defined as lim f(z) for z — > z% in the fundamental 
domain. 
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All the modes mediate a (gravitational) attraction. The sum ()D.20j) includes the 
long range contribution from the zero mode, giving rise to the standard Newtonian 
potential, characterized by the Planck mass, Mp ()80j) . while each Kaluza-Klein mode 
n provides a Yukawa- type force contribution, which is relevant at r ^ m~ l . We 
are interested in computing the gravitational potential at the "intermediate" dis- 
tances (|82p . For these distances, the largest contribution to the potential (|D.20|) is 
given by modes with mass m satisfying: 



A x k < m < k . 



(D.21) 



We can use this condition in the expansion of the Bessel functions which characterize 
the bulk profile of the Kaluza-Klein modes (|2U|) . Moreover, the mass spectrum in 
this range is given in Eq. ([67)1 . 

The "plus" and "minus" branches can be discussed simultaneously by introducing 
the parameters v± = ± (4a± — 1) . They satisfy 2 a/1 + a — v± =p 1 , so that v± ranges 
from ±1 to +oo , as a ranges from —1 to +oo (the RS point is at v_ — 1). We 
can set G\ = 1 in the mode solutions ()2()jl . since the normalization cancels in ()D.20j) 
(also, for shorthand we write v rather than v±). We then have: 



Z n (Z) = A- 



which give the following exact expressions: 
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(D.24) 



Inserting this ratio back into the potential (jD.20|) . and taking the expression (|8U|) 
for the 4D Planck mass on the UV brane, we obtain, 
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(D.25) 
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In the range (|D.21j) . the sum becomes: 
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where we have first used the result (|2()jl . and then we have approximated the expres- 
sion for the masses as m n ~ ukAin in the relevant regime. The potential (jD.25|) then 
rewrites: 

\ + kM 2 P 1 T{2\u ± \) 



T/ ^ 



M 2 r 



(D.27) 



M 3 (A;r) 2 l^±l 2 2 l^±l- 1 r(|z/ ± |) 2 _ ' 

where we have reintroduced the suffix ± which characterizes the two branches. 

The ratio kMp/M 3 controls the relative contribution of the zero mode and the 
Kaluza- Klein tower. For a < 1/4 , the zero mode is more localized towards the UV 
brane, and the massive Kaluza-Klein modes have a negligible effect in the regime (|82|). 
Substituting the value for M P given in eq. ()8L)j). and identifying z/_ = 1 — 4a = £ , we 
obtain the result (|83|) of the main text. In the complementary interval a > 1/4 , the 
zero mode is localized towards the IR brane, and the relative contribution from the 
Kaluza-Klein modes significantly increases. This interval is covered both by v_ < , 
and u + . We can describe it by a unique parameter £ = 4a — 1 , ranging from zero 
to infinity, which is identified with — 1/_ from to 1, and with v + from 1 to oo. This 
leads to the expression (|84j) given in the main text. 
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